Detection of Fermi Pairing via Electromagnetically Induced Transparency 
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An optical spectroscopic method based on the principle of electromagnetically-induced trans- 
parency (EIT) is proposed as quite a generic probing tool that provides valuable insights into the 
nature of Fermi paring in ultracold Fermi gases of two hyperfine states. This technique has the 
capability of allowing spectroscopic response to be determined in a nearly non-destructive manner 
and the whole spectrum may be obtained by scanning the probe laser frequency faster than the 
lifetime of the sample without re-preparing the atomic sample repeatedly. A quasiparticle picture is 
constructed to facilitate a simple physical explanation of the pairing signature in the EIT spectra. 

PACS numbers: 03.75.Ss, 05.30.Fk, 32.80.Qk 
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I. INTRODUCTION 

The two-component degenerate Fermi gas, in which 
the interaction between atoms of two different hyperfine 
states is made magnetically tunable via Feshbach reso- 
nance, has been the main source of inspiration for much 
recent excitement at the forefront of ultracold atomic 
physics research. In addition to being an ideal system 
for the exploration of the crossover from a Bose-Einsten 
condensate (BEC) of highly localized pairs to nonlocal 
Bardeen-Cooper-Schrieffer (BCS) pairs, the degenerate 
Fermi gas, when operating in the unitarity regime, con- 
stitutes a strongly interacting Fermi gas exhibiting a rich 
set of physics, the study of which may shed light on long- 
standing problems in many different branches of physics, 
in particular, condensed matter physics. 

A unique phenomenon of low temperature Femi system 
is the formation of correlated Fermi pairs. How to detect 
pair formation in an indisputable fashion has remained a 
central problem in the study of ultracold atomic physics. 
Unlike the BEC transition of bosons for which the phase 
transition is accompanied by an easily detectable dras- 
tic change in atomic density profile, the onset of pairing 
in Fermi gases does not result in measurable changes in 
fermion density. Early proposals sought the BCS pair- 
ing signature from the images of off-resonance scattering 
light The underlying idea is that to gain pairing in- 
formation, measurement must go beyond the first-order 
coherence, for example, to the density-density correla- 
tion. This is also the foundation for other detecting 
methods such as spatial noise correlations in the image 
of the expanding gas ^] , Bragg scattering [E 01 j Raman 
spectroscopy Q , Stokes scattering method [6| , radio fre- 
quency (RF) spectroscopy 't', 's'I , optical detection of ab- 
sorption 0, and interfcrometric method JXJJ. Among all 
these methods, RF spectroscopy 0, llj has been the only 
one implemented in current experiments [Til . [T^ . 

In this paper, we propose an alternative detection 
scheme, whose principle of operation is illustrated in 
Fig. [TJa). In our scheme, a relatively strong coupling 



and a weak probe laser field between the excited state 
e) and, respectively, the ground state \g) and the spin 
up state It), form a A-type energy diagram, which fa- 
cilitates the use of the principle of electromagnetically 
induced transparency (EIT) to determine the nature of 
pairing in the interacting Fermi gas of two hyperfine spin 
states: ||) and EIT [l^, in which a probe laser field 
experiences (virtually) no absorption but steep dispersion 
when operating around an atomic transition frequency, 
has been at the forefront of many exciting developments 
in the field of quantum optics [l4| . Such a phenomenon is 
based on quantum interference, which is absent in mea- 
surement schemes such as in Ref. @, where lasers are 
tuned far away from single-photon resonance. In the con- 
text of ultracold atoms, an important example is the ex- 
perimental demonstration of dramatic reduction of light 
speed in the EIT medium in the form of Bose conden- 
sate [iHl- This experiment has led to a renewed interest 
in EIT, motivated primarily at the prospect of the new 
possibilities that the slow speed and low intensity light 
may add to nonlinear optics [16] and quantum informa- 
tion processing 17:]. More recently, EIT has been used 
to spectroscopically probe ultracold Rydberg atoms [l^ . 
In this work, we will show how EIT can be exploited to 
reveal the nature of pairing in Fermi gases. 

Before we present our detailed calculation, let us first 
compare the proposed EIT method with the RF spec- 
troscopy method which is widely used in probing Fermi 
gases nowadays. In the latter 0, an atomic sample 
is prepared and an RF pulse is applied to the sample 
which couples one of the pairing states to a third atomic 
level 1 3). This is followed by a destructive measurement 
of the transferred atom numbers using absorption laser 
imaging. The RF signal is defined as the average rate 
change of the population in state |3) during the RF pulse, 
which can be inferred from the measured loss of atoms 
in I I). This process is repeated for another RF pulse 
with a different frequency. In addition to sparking many 
theoretical activities [l9|, US HH, [H, [H, this method 
has recently been expanded into the imbalanced Fermi 
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FIG. 1: (Color online) (a) The bare state picture of our model. 

(b) The dressed state picture of our model equivalent to (a). 

(c) A possible realization in ^Li. Here the states labelled by \i) 
{i = 1, 2, 6) are the 6 ground state hyperfine states. Most 
experiments involving ®Li are performed with a magnetic field 
strength tuned near a Feshbach resonance at 834G. Under 
such a magnetic field, the magnetic quantum number for the 
nuclear spin mi is, to a very good approximation, a good 
quantum number. The values of mi are shown in the level di- 
agrams. Two-photon transition can only occur between states 
with the same mi. Any pair of the lower manifold (|1), |2), 
and |3)) can be chosen to form the pairing states. In the ex- 
ample shown here, we choose |1) — jj), |2) = |J,) and |6) — \g). 
The excited state |e) (not shown) can be chosen properly as 
one of the electronic p state. 



gas systems [H, [2^ [l^, [H, [2^ , where paring can result 
in a number of interesting phenomena ^3Q] . A disad- 
vantage of this method is its inefficiency: The sample 
must be prepared repeatedly for each RF pulse. In addi- 
tion, for the most commonly used fermionic atom species, 
i.e., ^Li, the state |3) interacts strongly with the pairing 
states due to the fact that all three states involved has 
pairwise Feshbach resonances at relatively close magnetic 
field strength. This leads to so-called final state effect 
[23 | which greatly complicates the interpretation of the 
RF spectrum. 

In the FIT method, by contrast, one can directly 
measure the absorption or transmission spectrum of the 
probe light. Applying a frequency scan faster than the 
lifetime of the atomic sample to the weak probe field, the 
whole spectrum can be recorded continuously in a nearly 
non-destructive fashion to the atomic sample. Further- 
more, FIT signal results from quantum interference and 
is extremely sensitive to the two-photon resonance condi- 
tion. The width of the FIT transparency window can be 
controlled by the coupling laser intensity and be made 



narrower than Ep- As we will show below, this prop- 
erty can be exploited to detect the onset of pairing as 
the pairing interaction shifts and destroys the two-photon 
resonance condition. In addition, due to different selec- 
tion rules compared with the RF method, one can pick 
a different final state whose interaction with the pairing 
states are negligible [see Fig. lUc)], hence avoiding the 
final state effects. 

The paper is organized as follows. In Sec. HIl we de- 
scribed the model under study and define the key quan- 
tity of the proposal — the absorption coefficient of the 
probe light. In Sec. IIIIl we present the expression of the 
probe absorption coefficient and construct a quasiparti- 
cle picture that will become convenient to explain the 
features of the spectrum. The results are presented in 
Sec. IIVI where spectral features at different temperatues 
are explained. We also show that how FIT spectrum can 
be used to detect the onset of pairing. A breif summary 
is presented in Sec. [V] Finally, we provide an appendix 
in Sec. IVII where the derivation of the FIT spectrum is 
provided. In particular, we include in this derivation the 
pairing fluctuations in the framework of the pseudogap 
theory jT|. 

II. MODEL 

Let us now describe our model in more detail, begin- 
ning with the definition of LUi and f2i as the temporal and 
Rabi frequencies of the probe (i = p) and coupling (i = c) 
laser field of plane waves copropagating with an almost 
identical wavevector (along z direction). The system 
to be considered is a homogeneous one with a total vol- 
ume V, and can thus be described by operators a\^,i (aj^ ■) 
for annihilating (creating) a fermionic atom in state \i) 
with momentum /ik, and kinetic energy ek — lT?k^/2m, 
where m is the atomic mass. Here, (ik.i are defined in 
an interaction picture in which Ok.e ~ oj,. ^e"*'^''*, ak,g 

= aL,3e'("<=""''')*, and du,a = =t, i), where aj^ - 

are the corresponding Schrodinger picture operators. 

In a probe spectrum, the signal to be measured is the 
probe laser field, which is modified by a polarization hav- 
ing the same mathematical form as the probe field ac- 
cording to [31] 

dz ^ c dt 2 ' 

where Pp is the slowly varying amplitude of that polar- 
ization, dij is the matrix element of the dipole moment 
operator between states \i) and and fiQ and c are 
the magnetic permeability and the speed of light in vac- 
uum, respectively. The parameter a in Eq. Ilj represents 
the complex absorption coefficient of the probe light 3l|. 
By performing an ensemble average of atomic dipole mo- 
ment, we can express a as 

" = ^?^E(4A+''-^>^'^'"''^ (2) 

P k,q 
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where ao = iJ,oLijpc\dei\ ■ The real and imaginary part 
of a correspond to the probe absorption and dispersion 
spectrum, respectively. 

To determine the probe spectrum, we start 
from the grand canonical Hamiltonian H — 

Z^k ('^ik + ■^2k + Hsk) , where 

Wlk = (e'fe - (5p) ak.e«k,e + (4 - 5) 4,gak,g , 

U2V = -^(^^cak+kz,,e«k.3 + f^pol+kj, ,e«k, T ) - h.C , 

^3k XI 4alc,o-«k,o- - (Aaj, -fat_i^^^ + h.c) , 

describe the bare atomic energies of states |e) and \g), 
the dipole interaction between atoms and laser fields, 
and the mean-field Hamiltonian for the spin up and 
down subsystem, respectively. Here, ej. = e/j — /l( with 
/i being the chemical potential, 5p = h{u}p — uje]) and 
5c ~ h{ujc ~ uJeg) Siie the single-photon detunings, and 
i5 = Sp — 6c is the two-photon detuning with coij be- 
ing the atomic transition frequency from level \i) to \ j). 
In arriving at Tisk, in order for the main physics to be 
most easily identified, we have expressed the collisions 
between atoms of opposite spins in terms of the gap pa- 
rameter A = ~UV~^ J^kiO'-'kd'^'k.A) under the assump- 
tion of BCS paring, where U characterizes the interaction 
between | f) and | j) which, in the calculation, will be 
replaced in favor of the s-wave scattering length as via 
the regularization procedure: 

m 11 X ^ 1 

4^t¥^s ^ U ^ V ^2^' 

k 

A more complex model including the pseudo-gap physics 



[19| will be presented later in the paper. Finally, we note 
that the effect of the collisions involving the final state 
|g) in the RF spectrum has been a topic of much recent 
discussion [2ll l22l In our model, the spectra are 

not limited to the RF regime, and this may provide us 
with more freedom to choose \g) (and |e)) that minimizes 
the final state effect. In what follows, for the sake of 
simplicity, we ignore the collisions involving states \g) 
(and |e)). In practice, the effects of final state interaction 
can be minimized by choosing the proper atomic species 
[s^ l or hyperfine spin states [33] . In the example shown in 
Fig- [He), it is indeed expected that \g) does not interact 
strongly with either of the pairing state. 

III. QUASIPARTICLE PICTURE 

The part of the Hamiltonian describing the pairing of 
the fermions can be diagonalized using the standard Bo- 
goliubov transformation: 

ak,T = Wfcdk.T + ^fe"Lk,i ' 

where Uk = ^{Ek + e'^) /2Ek,Vk = ^ {Ek - e'^) /2Ek, 
and Ek = e'^ + A'^ is the quasiparticle energy disper- 
sion. Now we introduce two sets of quasiparticle states 
± Ik), representing the electron and hole branches, re- 
spectively. The corresponding field operators are defined 
as 

Q!k,+i = dk,T , Q!k,-i = aLkj J 

in terms of which, the grand canonical Hamiltonian can 
be written as 



H = Y1 [^'^'k ^ ^p) «k,eak,e + (4 - S) al gUk^g + £^fcak,+iak,+i - £^/cak -i^k. 



2 "k+kL,e"k,g 



h.C. - 



-^ak+ki,,e"k.-l 



(3) 



A physical picture emerges from this Hamiltonian very 
nicely. The state |-l-lk) (|"lk)) has an energy dispersion 
+Ek (—Ek) and is coupled to the excited state |e) by 
an effective Rabi frequency flpUk {flpVk), which is now 
a function of k. In the quasiparticle picture, our model 
becomes a double A system as illustrated in Fig. [Ijb). 
Let -|-A (—A) denote the A configuration involving |+lk} 
(| — Ik))- The -fA (—A) system is characterized with a 
single-photon detuning of Sp + Ek {Sp — Ek) and a two- 
photon detuning of S + E^ {S — Ek). In thermal equi- 
librium at temperature T (in the absence of the probe 



field), we have 

(al _,_iak'.+i) = (^k.k' - {al _iak'.-i) = ^k.k'/ {Ek) , 

(4) 

where 

f{Lu)^[e^p{uj/kBT) + l]-\ (5) 

is the standard Fermi-Dirac distribution for quasiparti- 
cles. Thus, as temperature increases from zero, the prob- 
ability of finding a quasiparticle in state |+lk) increases 
while that in state | — Ik) decreases but the total proba- 
bility within each momentum group remains unchanged. 
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Similarly, in the quasiparticle picture, the probe spec- 
trum receives contributions from two transitions 

k.q 

[UqPe, + l (k + ki, q) + VqPe-1 (k + k^, q)] , (6) 

where pi^±i(k, k') = (^(xl^i ±ia]i,ij are the off-diagonal 
density matrix elements in momentum space. 

The equations for the density matrix elements can be 
obtained by averaging, with respect to the thermal equi- 
librium defined in Eq. the corresponding Hcisen- 
berg's equations of motion based upon Hamiltonian (3). 
In the regime where the linear response theory holds, the 
terms at the second order and higher can be ignored, and 
the density matrix elements correct up to the first order 
in rip are then found to be governed by the following 
coupled equations: 



in — 
dt 



Pe^^ (k + kL,q) 
Pg-n (k, q) 



= M„ 



/5e,.^(k + ki,q) 
Pg,v (k, q) 



- (k) <5k.q, iv 



±1) 



(7) 
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FIG. 2; (Color online) (a) A (black solid curve) and the probe 
absorption coefficient real(a) at 5 = (red dotted curve) as 
functions of T, obtained from the mean-field BCS theory, (b) 
Real(a) as a function of S (absorption spectrum) at different 
T. (c) A, Asc and Apg as functions of T obtained from the 
pseudogap approach. (Asc ~ and A = Apg when Tc < T < 
T*). (d) A comparison of various absorption spectra at T = 
O.STf. The parameters are 5^ = 0, 7 = 3S0Ef (~ lOMHz), 
Q.C = 5Ef (~ O.lMHz), and l/ikpas) = -0.1. 
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Here we have introduced phenomcnologically the pa- 
rameter 7 which represents the decay rate of the ex- 
cited state |e). Inserting the steady-state solution from 
Eq. ([7]) into Eq. ([S]), we immediately arrive at a (Jc, S) = 
a+i ((5c, <5) + q;_i ((5c, <5), where 



±i(,5c,(5) = *^5]«;k(^c,'5,±ii;fc)/(±£;;c) j g , 



with 

Wk (^c,<5,a;) 



Ak (^c,<5,a;) (e'^ - (5 - c^) - | % 



(8) 

71^ (9) 



and Ak ((5c, (5, w) = £^+^1. 



■ (5 — 17 — w. 



IV. RESULTS 

Examples of the probe absorption coefScient, Re(a), 
are presented in Fig. [DJa) and (b). For the results 
shown in this paper, we choose l/(fcj7as) = —0.1 where 
we denote Ep^kp, and Tp = Ep/ks be Fermi energy. 



wavenumber, and temperature, respectively, for the non- 
interacting Fermi gas. The black solid line in Fig. [2ta) 
represents the gap parameter in the mean-field calcula- 
tion, from which we can see that the critical tempera- 
ture below which the system exhibits pairing is given by 
Tc — 0.435Tp for the parameters chosen. The dotted red 
curve in Fig.[2{a) represents the absorption coefficient at 
two-photon resonance (5 = 0. We can see that it remains 
at zero for T > Tc but increases sharply once the tem- 
perature drops below Tc. We note that this feature can 
be used as a sensitive gauge for detecting the onset of 
Fermi pairing. With this being emphasized, we now turn 
to explain the main spectroscopic features displayed in 
Fig.EIb). 

First, as long as T > Tc where A = 0, one can show 
that the spectrum is essentially independent of T and 

Re(Q;) oc 



[(,5 + (5c)(5 - |f^c/2|2]2 + ,5272 



From this expression, one can easily see that there 
exists around (5 = a narrow transparency window 
with a width determined by the optical pumping rate 
Fop = 1 ric I ^ 7/ [4 ((52 -I- 72)] [see the blue dashed curve 
for T = 0.5Tp in Fig. El^b)]. This feature can be most 
easily understood from the bare state picture [Fig. [Ha)], 
where state |t) is decoupled from state |J,) so that the 
spectrum is of FIT type for a A system involving |e), 
\g), and |t). Further, because states \g) and |t) share 
the same energy dispersion ej, , the two-photon resonance 
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condition (5 = holds for atoms of any velocity groups; 
the absence of absorption at i5 = signals the existence 
of a coherent population trapping state. 

As T decreases below Tc, a double-peak structure de- 
velops [see the red dotted hue for T = OATp in Fig.IlJb)]. 
The two peaks can be understood as contributed by the 
quasiparticle state | -I- Ik) and | — Ik), respectively. In 
the limit where T is far below Tc [see the black solid 
line for T = O.OITf in Fig. [^b)], +A system has negli- 
gible contribution to the probe spectrum because there 
exists virtually no quasiparticles in state [-l-lk)- Thus, 
the spectrum is solely contributed by —A system, result- 
ing a single-peak structure. However, unlike the situa- 
tions above Tc, here while the dispersion of an atom in 
state \g) continues to be ej,., the dispersion of a dressed 
particle in state | — Ik) is —Ek- As a result, the effec- 
tive two-photon resonance condition e'^, — 6 + — 
is now momentum dependent. Aside from a shift, the 
transparency window becomes inhomogeneously broad- 
ened with a linewidth in the order of Ep- A consequence 
of the momentum-dependence of the two-photon reso- 
nance condition is that, for any given probe laser fre- 
quency, only atoms with the 'right' momentum result in 
perfect destructive quantum interference. Consequently, 
Re(a) can no longer be zero for any probe frequency. 
This underlies the sharp increase of the probe absorption 
at 5 = below Tc as shown in Fig. [D^a). 

We also want to emphasize that the spectrum shown 
in Fig. Eljb) can be obtained by scanning the probe laser 
frequency over a range on the order of Ep ~ O.lMHz. 
We may take typical spectral features of the Fermi gas 
to be Su) ~ Q.lEp ~ lOKHz. To resolve such features, 
using the energy-time uncertainty relation, we can use 
a scan rate of lOKHz/O.lms, then the total scan time 
can be estimated to be around 1 ms. As this time is 
much shorter compared with the typical lifetime of the 
Fermi gas, this method can be regarded as nearly non- 
destructive. This demonstrates the great efficiency of the 
KIT probe. 

In a more realistic model where pair fluctuations are 
included, gap A is divided into a BCS gap A^c for con- 
densed (BCS) pairs below Tc and a pseudogap Apg for 
preformed (finite momentum) pairs below temperature 
T* according to A^ = A^^ + ^Ig O- Results includ- 
ing pseudogap physics are illustrated in Fig. [DJc) and 
(d) and the detailed derivations can be found in the Ap- 
pendex. In contrast to the weakly interacting regime, 
where T* is virtually the same as Tc, T* is much higher 
than Tc in strongly interacting regime as is clearly the 
case of present study according to Fig. [H^c). It needs to 
be stressed that pair fluctuations can result in a finite 
lifetime 7"^ for preformed pairs which tend to broaden 
the spectral features, so that only when 7p is sufficiently 
small can the double-peak spectroscopic structure be re- 
solved as Fig.[2l[d) demonstrates. Finally, the two-photon 
resonance here is only sensitive to A because Ek de- 
pends on the total gap A [l^. As a result, like its RF 
counterpart [1^, the KIT method cannot distinguish be- 



tween Asc and Apg. However, the qualitative features 
of Fig. are not changed as long as we regard the 
corresponding critical temperature as T* . 



V. SUMMARY 

In summary, we propose to use optical spectroscopy in 
an FIT setting to probe the fermionic pairing in Fermi 
gases. We have demonstrated that the FIT technique of- 
fers an extremely efficient probing method and is capable 
of detecting the onset of pair formation (i.e., determin- 
ing T*) due to its spectral sensitivity. With a sufficiently 
weak probe field, the whole spectrum may be obtained 
with a nearly non-destructive fashion via a relatively fast 
scan of probe frequency, without the need of repeatedly 
re-preparing the sample. We note that in this work, we 
have focused on probing the atomic system using pho- 
tons. In the future, it will also be interesting to study 
how we can use atomic Fermi gas to manipulate the light. 
Superfluid fermions can serve a new type of nonlinear me- 
dia for photons. Finally, we want to remark that, in this 
work, as a proof-of-principle, we have only considered a 
homogeneous system. As usual, the trap inhomogeneity 
can be easily accounted for within local density approxi- 
mation. Nevertheless, we note that the capability of de- 
tecting the onset of pairing remains the same even in the 
presence of the trap. Furthermore, as optical fields are 
used in this scheme, one may focus the probe laser beam 
such that only a small localized portion of the atomic 
cloud is probed, hence there is no need to average over 
the whole cloud. 
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VI. APPENDIX: EIT SPECTRA INCLUDING 
PSEUDOGAP 

In this appendix, we generalize the result of Eq. ([2]) 
for a valid under the mean-field BCS pairing to a more 
realistic situation where pair fiuctuations are included in 
the form of pseudogap. We show two different ways to 
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accomplish this generahzation. The first is an approach 
used more often by people working in the field of quantum 
optics. The second uses the linear response theory Q 
more familiar in the field of condensed matter physics. 



A. A Brief Account of Pseudogap Theory 

First, let us highlight the results of pseudogap theory 
[l9| that are relevant to our EIT spectrum calculation. 
When pairing fluctuations at finite temperature are in- 
cluded in the framework of the pseudogap model [l^ , the 
BCS gap equation and number equation are still valid. 
However, the gap A is now regarded as the total gap 
divided into a BCS gap Age for condensed (BCS) pairs 
below Tc and a pseudogap Apg for preformed (finite mo- 
mentum) pairs: 



= Al. 



Ai 



The onset of the total gap A occurs at temperature T* , 
which is greater than Tc. The system with preformed 
pairs is described by the Green's function 

G-i(k,zm„) = G^i(k,iu;„) - S(k,iu;„), (10) 

where the non-interacting Green's function 

Gq ^(k,iu;„) = (iw„ - ei^)~\ (11) 

and the self energy 

E(k, iwn) = Ssc(k, iwn) + Spg(k, iw^) 

A^ A^ 
= ""-h— ^ (12) 

%Wn + Cfe + f/c + *7p 

with Wn being the fermi Matsubara frequency and 7"^ 
the finite lifetime of pseudogap pairs. The spectral func- 
tion A(k, w) can be obtained from the Green's function 
via the relation 

A(k, cj) = -2 Im G (k, + iO+) , 

which, with the help of Eqs. ([TOl), (HIl), and ([T2|), is found 
to be given by 



2(u; + .;)^7pAg^ 



(13) 



where £:f = A/e',. + AL. In the limit of 



7p 



and 



^ Ek,we recover from Eq. (fT3|) the spectral function 
under the BCS paring 

A{k, w) = 2tt[uIS{lj - Ek) + vl5{u + Ek)] . (14) 



terms of the equal time correlation function /iq.k (i) 
^ (i)ak+ki,,e {t)) as 



"q.T 



.ao 1 ,. 

\lp V t ►oc 



k,q 



{t)e 



i(k~q) -r 



(15) 



where the limit is introduced to indicate explicitly that 
we are interested in the steady state spectrum. Here, 
dq I (t) and dk+k^.e (^) obey the Heisenberg equations of 
motion 



ih 



dt 



with 



flk+ki, 



M 



(16) 



(17) 



Note that due to the dissipative nature of our model, 
strictly speaking, Eqs. (fTB|) should be those of quantum 
Langevin equations containing the noise operators of the 
reservoir that gives rise to the decay rate 7. Here, in 
anticipation that Eqs. (fTB]) will produce the right aver- 
ages of our interest, we have ignored the noise operators. 
We solve Eqs. (IT51) for dk+kj,.e (t) in the limit of t ^ cxo 
when the terms involving the initial operators have all 
died away, and then combine it with d^ ^ (t) to form 



it) = ^ 



-jA/(t-t') 



G< (k,i',i)(5k,qdt'. 



(18) 

where denotes the element at the first row and the 

first column of the matrix inside the square bracket, and 

G< (k, <',i) = i ^dj^ ^ (t) 0-^,1 (^')^ is one of the Green's 
functions in real time. By substituting G^ (k, t' , t) in 
Eq. (jlSp with a Fourier transformation of its counterpart 
in real frequency, G^ (k, uj), we are able to carry out the 
time integration in Eq. p8|) explicitly, leading to 



3k, q- 



A(k,u;)fiLo) 

M ~LU 



where the use of a well-known relation: G^ (k, to) = 
if (w) A (k, Lj) [s^ has been made. Finally, replacing 
[1/{M — a;)]ii with Wk i^ci^,^): obtained with the help 
of Eq. pT]) , we arrive at 



" = *^E/ ^A{Kc.)f{cj)w^{S,,S,u;) (19) 



B. Quantum Optics Approach 

In order to develop a formalism which directly incor- 
porates the spectral function, we rewrite Eq. ([5]) in 



where Wk (Sc, uj) is defined in Eq. ([9]) of the main text. 
One can easily check that Eq. (fT9|) reduces to Eq. ([8]) in 
the limit of mean-field BCS pairing when Eq. (fTi|) is used 
as the spectral function. 
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C. Condensed Matter Approach 

In order to use the linear response theory widely used 
in condensed matter physics, we first divide our sys- 
tem into a "left part" comprising two hyperfine spin 
states: |t) and ||), whose physics has been described in 
Sec. IVI A[ a "right part" consisting of the coupling laser 
field and states \g) and |e), described by the Hamiltonian 



Hn = 



EN 



2 E^k+ki.e^M + '^-C- 

k 



and finally the coupling between the two parts induced by 
the probe field, described by the tunneling Hamiltonian 

= -^E^k+k^.e^kj + h.c 

k 

= A + Al 

Next, we change Hr into a diagonal form 



(20) 



in terms of a pair of dressed state operators, a and (3, 
defined via the transformation 



a'k+ki,,e 
«k,g 



,/3 



/3k 



(21) 



where 



^ [Ck + ??fe ± \/(Cfc 



rik) 



(22) 
(23) 



with C,k = ^'k+kL ~ ^'^^ ~ ^'k ~ ^- terms of the 
dressed state operators, A becomes 



(24) 



and is in a form to which the linear response theory (3j| 
is directly applicable. Following the standard practice, 
we then find 



1_ O J 



k r]—a,f3 



+ 00 



2n 



(25) 



In Eq. ((25| . j4i(k, cj^) is same as A(k, w^) defined in 
Eq. din]), while A'^(k, ujr) is given by 2n6 {ujr - E^) be- 
cause the right part is in a normal state described by the 
Green's function (k, iw„) — iwn — E^^. Integrating 
over lor, we change Eq. (j25p into 



k rj—OL^P 



+ 00 

2 / 



27r 



^(k,w 



(26) 

where the dummy variable ljl has been changed into lj. 
We now include the effect of the decay of the excited 
state phenomenologically by replacing 5p with 5p — i^y. 
We see that E^ now become imaginary which signals the 
inability of the dressed states to hold populations. This 
along with the fact that the dressed states here are the 
superpositions of the initially empty states provide us 
with the justification to set /(i^^) = in Eq. (QH). With 
these considerations, we finally arrive at 



= -T E / u;)fiu;)w^ (<5„ S, u;) (27) 

k „ 



where the use of Eqs. ((22|) and ((23)) is made. It is clear 
from Eq. (fT9|) that a is proportional to i{A) in Eq. ([27]) . 
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